A search for light dark matter in XENON10 data 
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We report results of a search for light (< 10 GeV) particle dark matter with the XENON10 
detector. The event trigger was sensitive to a single electron, with the analysis threshold of 5 
electrons corresponding to 1.4 keV nuclear recoil energy. Considering spin-independent dark matter- 
nucleon scattering, we exclude cross sections a n > 7 x 10~ 42 cm 2 , for a dark matter particle mass 
m x = 7 GeV. We find that our data strongly constrain recent elastic dark matter interpretations 
of excess low-energy events observed by CoGeNT and CRESST-II, as well as the DAMA annual 
modulation signal. 



Recently, the CRESST-II and CoGeNT collaborations 
have reported the observation of low-energy events in 
excess of known backgrounds [TJ [2]. This has encour- 
aged the hypothesis that these signals — in addition to 
the long-standing DAMA [3] annual modulation signal 
— might arise from the scattering of a light (< 10 GeV) 
dark matter particle |4UT0]. The CDMS [HI H] dark 
matter search data been have re-analyzed with lowered 
energy thresholds, but do not fully exclude light dark 
matter interpretations. In order to maintain sensitiv- 
ity to the scattering of such light galactic dark matter, 
an experiment needs either a very low O(keV) energy 
threshold, as with CoGeNT, or light target nuclei, such 
as the oxygen atoms in the CRESST-II detector. This 
is a simple requirement of the kinematics [13] . For ex- 
ample, consider a halo-bound dark matter particle with 
mass m x = 10 GeV and velocity 600 km s _1 . The maxi- 
mum recoil energy that would result from an elastic scat- 
ter of such a particle in an earth-bound target would be 
about 20 keV for a recoiling oxygen nucleus (CRESST- 
II), 9 keV for germanium (CoGeNT) and only 6 keV for 
xenon. The respective energy thresholds of CRESST-II 
and CoGeNT are approximately 10 keV [T3] and 2 keV 
[5] , which combined with a low background event rate re- 
sults in good sensitivity to light mass dark matter. The 
energy threshold of previously reported XENON10 data 



[151 116] depended on the primary scintillation efficiency 
of liquid xenon for nuclear recoils (£ e //) (TTJ [TB]. For 
a conservative assumption of the energy dependence of 
£eff [H]j the threshold was about 5 keV. 

It is possible to obtain a lower energy threshold from 
existing XENON10 dark matter search data, if nuclear 
recoil energy is measured by the detected electron sig- 
nal. The method allows us to reach an energy threshold 
E„ r ~ 1 keV. At such low nuclear recoil energies, the pri- 
mary scintillation signal is generally absent. As a result, 
two important aspects of the XENON10 detector [TTJ] 
performance are compromised: the ability to precisely 
reconstruct the z coordinate of a particle interaction, and 
the discrimination between incident particle types. The 
detected ratio of scintillation to electron signals was used 
in jTS], fTB] to discriminate and reject about 99.5% of elec- 
tromagnetic background events which would otherwise 
have been treated as dark matter candidate events. Loss 
of this discrimination thus reduces our sensitivity to nu- 
clear recoils from dark matter scattering, by about two 
orders of magnitude. Still, the lower energy threshold 
we obtain permits the exploration of new regions of dark 
matter particle mass (m x ) and cross section (a n ) param- 
eter space. 

The XENON10 detector is described in detail in [15] . 
For the present discussion, we remind the reader that 
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XENON10 observes particle interactions via detected pri- 
mary scintillation photons (the SI signal) and electrons 
(the S2 signal). The electrons are drifted across the liq- 
uid xenon target and extracted into gaseous xenon, where 
they create proportional scintillation photons. This al- 
lows both SI and S2 to be observed with the same photo- 
multiplier tube arrays. A single electron extracted from 
the liquid target results in about 27 photoelectrons in the 
photomultiplier arrays. It is this robust signal that gives 
the S2 channel its lower energy threshold. The (x, y) co- 
ordinates of particle interactions are reconstructed from 
the hit pattern of the S2 signal on the top photomulti- 
pliers. The z coordinate is usually reconstructed from 
z = VdAt, where Vd — 0.20 cm /zs _1 is the electron drift 
velocity [23] and At is the measured time delay between 
the SI and S2. This method cannot be used if there is 
no SI signal, as is often the case for very low energy nu- 
clear recoils [24]. The S2 pulse width a e carries a mild 
^-dependence [25] . Larger detectors [26l [27] , if operated 
with a lower value of E^, may obtain a reliable estimate 
of the z coordinate from <r e [2B]; this should lead to a 
significant rejection of edge (in z) events and a commen- 
surate improvement in sensitivity. 

Ideally we would like to reconstruct the nuclear recoil 
energy for each event from E„ r = e(n 7 +n e )/f n , as in 
[53] [3D], with e = 13.8 eV the average energy to create a 
photon or electron, and /„ the nuclear recoil quenching 
[31 j . However, at low recoil energies the small number n 7 
of primary scintillation photons often does not result in a 
measurable SI response. Since we are interested in events 
at very low recoil energy, we calibrate the energy scale 
using only the number n e of measured electrons in the 
S2 signal. The electron yield of liquid xenon for nuclear 
recoils has been measured directly using tagged neutron 
scattering [T7| 132] . The lowest-energy data point from 
[T7] implies an S2 signal of 36±6 electrons in the 3—5 keV 
range, as shown in Fig. [T] An indirect measurement of 
the electron yield is described in 25 , and was obtained 
following the method detailed in [21]. The central and 
±ler contours of that work are shown in Fig. [I] (dash-dot 
curves). That Q y rises with decreasing E„ r between 100 
and 10 keV is a result of the increasing fraction of nuclear 
recoil energy given to electrons (rather than photons) 
over that energy range [21] ■ 

Following [29] , we obtain a theoretical prediction for 
the electron yield, 

In this equation, N ex /Ni is the number ratio of excited 
to ionized xenon atoms, and £ = Nia/(4a 2 v) is the sin- 
gle parameter upon which the Thomas-Imel box model 
depends j33j [34]. In Eq. [I] we explicitly indicate the 
dependence of f n on the proportionality constant k, be- 
tween the velocity of a xenon nucleus and its electronic 
stopping power. In Fig. [I] we show Eq. [T] predictions 
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FIG. 1. The electron yield Q y of liquid xenon for nuclear re- 
coils. Theoretical curves (solid and dashed) were calculated 
from Eq. [T] as described in the text. Also shown are measure- 
ments from [17] (~k), [25] (dash-dot curve, with ±lcr contours) 
and [32] (O and ■, uncertainty omitted for clarity). 



for two k values, which correspond to calculations by 
Lindhard [3T] and Hitachi [32]. The dashed curve is the 
best-fit case from j2FJ, with k = 0.166, N ex /Ni = 1.05 
and 4£/iVi = 0.024. The solid curve, which we will use to 
calibrate the energy scale in the present work, takes the 
more conservative k = 0.110, from which we obtain the 
best-fit parameters N ex /Ni = 1.09 and 4£/A t = 0.032. 
This results in the most conservative a n exclusion limits 
based on available data and theoretical considerations, 
and is consistent with our neutron calibration data [25) . 
However, it is in tension with the measurements of Ref. 
[17] below ~ 7 keV. As discussed in [35], the rising mea- 
sured Q y values in this regime could be influenced by 
trigger threshold bias. We emphasize that our energy 
calibration (Fig. [I] solid curve) relies on extrapolation of 
Lindhard's theory [3T] below 4 keV. 

We report results from a 12.5 live day dark matter 
search, obtained between August 23 and September 14, 
2006. This data set is distinct from the previously re- 
ported PS] [H] XENON10 dark matter search data: the 
present data was obtained with the the secondary scintil- 
lation gain about 12% higher, and the S2-sensitive trigger 
threshold set at the level of a single electron. The trigger 
efficiency for single electrons is > 0.80 [56] . 

Event selection criteria are summarized in Table [I] 
Candidate events were required to have an S2 signal of 
at least 5 electrons, or 1.4 keV. Although the detector 
is sensitive to smaller energy depositions, we are unable 
to reliably assess the cut acceptance e c for smaller nu- 
clear recoils, due to the trigger configuration during the 
neutron calibration 19 . The position of interaction was 
required to fall within r < 3 cm. This central region 
features optimal self-shielding by the surrounding xenon 
target. A signal-to- noise cut required the S2 pulse to 
contain at least 0.45 of the total area of the event record. 
The acceptance of this cut rises monotonically from 0.94 
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to > 0.99 between 1.4 keV and 10 keV. Valid single scat- 
ter event records were required to have only a single S2 
pulse of size > 4 electrons. Events in which an SI signal 
was found were required to have logio(S2/Sl) within the 
±3(7 band for elastic single scatter nuclear recoils. This 
band was determined from the neutron calibration data, 
and has been reported in a previous article |15j . Events 
in which no SI signal was found were assumed to be dark 
matter candidate events and were retained. 



TABLE I. Summary of cuts applied to 15 kg-days of dark 
matter search data, corresponding acceptance for nuclear re- 
coils e c and number of events remaining in the range 1.4 < 
E nr < 10 keV. 





Cut description 


e c 




1. 


event localization r < 3 cm 


1.00 a 


125 


2. 


signal-to-noise 


> 0.94 


58 


3. 


single scatter (single S2) 


> 0.99 


38 


4. 


±3cr nuclear recoil band 


> 0.99 


23 



a limits effective target mass to 1.2 kg 

The remaining events in the lowest-energy region are 
shown in Fig. [2] (left) versus their S2 pulse width a e . The 
equivalent number of electrons is indicated by the inset 
scale. Events in which an SI signal was observed are indi- 
cated by a circle. Figure [2] (right) shows the width profile 
of the S2 signal in the top, middle and bottom third of 
the detector, based on single scatter nuclear recoils with 
known At and 5 < S2 < 100 electrons. Gaussian fits are 
shown to guide the eye. 

The top panel of Fig. [2] shows the distribution of re- 
maining candidate events (+) with S2 < 4 electrons. The 
distribution of background single electron events, sam- 
pled from a time window at least 20 /is after higher- 
energy events, is also shown (A). The single electron 
background events are a subject of ongoing study, and ap- 
pear to originate from multiple physical phenomena. One 
possibility involves photoionization of impurities in the 
liquid xenon |37j . Another possible origin is from excess 
free electrons trapped at the liquid surface. This could 
occur because the emission of electrons from the liquid 
to the gas is nearly — but likely not exactly — unity [35] . 
As a result, every S2 signal could be a potential source of 
a small number of trapped electrons. Delayed emission 
of the trapped electrons may result from the requirement 
that both the electron kinetic energy and the z compo- 
nent of the electron momentum be sufficient to overcome 
the surface potential barrier [3"9"] . 

The signal-to-noise cut was motivated by a distinct but 
closely related class of background event, which consists 
of a train of approximately ten to several tens of single 
electrons over a period of O(100 fis). The origin of these 
events is also not yet clear. Often several single electrons 
in an electron train overlap in time, to the degree that 
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FIG. 2. (left) All candidate dark matter events remaining 
( x and O) after the cuts listed in Table I] Events in which an 

51 was found are shown as O- The number of electrons in the 

52 signal is indicated by the inset scale, (top) Distribution 
of candidate events with < 4 electrons (+), and distribution 
of background single electrons (A) as described in the text, 
(right) S2 pulse width distributions for single scatter nuclear 
recoils in the top, middle and bottom third of the detector. 



they appear as a single S2 pulse containing ^2 — 6 elec- 
trons. These spurious pulses often have a e > 0.30 (the 
3a width for a single electron) and so could be removed 
based on pulse width. However, the signal-to-noise cut 
more precisely targets the presence of multiple additional 
single electrons in the event record. 

The energy resolution for S2 signals depends primarily 
on Poisson fluctuation in the number of detected elec- 
trons, with an additional component due to instrumen- 
tal fluctuations. This is discussed in detail in [35) . and 
for higher energy signals in [19]. So as not to over- 
state the energy resolution, we adopt a parameteriza- 
tion which follows the Poisson component only, given by 
lZ(E nr ) = (2E nr ) -1 / 2 . We assume a sharp cutoff in Q y at 
E.,„. = 1.4 keV, and then convolve the resolution with the 
predicted differential dark matter scattering rate. This 
ensures that a n exclusion limits are not influenced by 
lower-energy extrapolation of the detector response. The 
scattering rate as a function of nuclear recoil energy was 
calculated in the usual manner [Q] (c/. [13). We take 
the rotational speed of the local standard of rest and 
the velocity dispersion of the dark matter halo to be 
Vq = 230 km s _1 , and the galactic escape velocity to be 
v esc — 600 km s _1 [H]. We use the p ma x method [32] to 
calculate 90% C.L. exclusion limits on the cross section 
a n for elastic spin-independent dark matter — nucleon 
scattering as a function of m x . All remaining events in 
the the range E„ r > 1.4 keV are treated as potential dark 
matter signal. The results are shown in Fig. [3j If Q y 
were 40% higher (lower) below 4 keV, the exclusion limits 
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would be about x2 stronger (weaker) at m x — 7 GeV. 
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FIG. 3. Curves indicate 90% C.L. exclusion limits on spin- 
independent o n for elastic dark matter scattering, obtained by 
CDMS (dotted [11], and dashed [12]) and XENON100 (dash- 
dot [5S]). The region consistent with assumption of a positive 
detection by CoGeNT is shown (hatched) [2], and (shaded) 
0J; the latter assumes a 30% exponential background. Also 
shown is the 3a allowed region for the DAMA annual modu- 
lation signal (solid contour) [40] . 

The exclusion limits and allowed regions shown in Fig. 
[3] assume a simple Maxwell-Boltzman distribution for the 
dark matter halo. Given the likelihood of significant de- 
partures from this distribution [43 , it is important to 
understand if astrophysical uncertainties could alter the 
incompatibility of our results with the positive detec- 
tion scenarios shown in Fig. [3] A method for doing 
so is described in [44 j . and predicts that not less than 
~ 5 counts keV" 1 kg -1 day" 1 should be observed in 
a xenon detector, if the unexplained low-energy rise ob- 
served by the CoGeNT detector [5] were due to dark mat- 
ter scattering. It can be seen from Table[l]that we observe 
an event rate of ~ 0.2 counts keV -1 kg -1 day -1 on the 
interval 1.4 < E m . < 10 keV, indicating that the order 
of magnitude exclusion of the CoGeNT regions shown in 
Fig. [3] is robust against astrophysical uncertainties. Due 
to the preliminary nature of the CRESST-II results we 
do not show a corresponding allowed region, although it 
appears likely to lie above the DAMA region, as shown 
in Fig. [3] of Ref. @5]. 

We have shown for the first time that it is possible to 
perform a sensitive search for dark matter with a liquid 
xenon time-projection chamber, using only the electron 
signal. The advantage of this analysis is an increased 
sensitivity to light (< 10 GeV) dark matter candidate 
particles, due to the approximate factor x5 decrease in 
the detector energy threshold. For larger particle masses, 
standard analyses [THl UHl EE] offer superior sensitivity. 
The present work appears to severely constrain recent 
light elastic dark matter interpretations of the excess low- 
energy events observed by CoGeNT and CRESST-II, as 
well as interpretations of the DAMA modulation signal. 
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